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Aims. We investigate the dependence of galaxy clustering on luminosity and stellar mass, in the redshift range 0.5 < z < 1.1, using the first 
~ 55000 redshifts from the VIMOS Public Extragalactic Redshift Survey (VIPERS). 
r ) • Methods. To characterize the spatial properties of the VIPERS galaxies, we measure the redshift-space two-point correlation functions (2PCF), 
' and ^{rp,n) , and the projected correlation function, Wp(rp), in samples covering different ranges of B-band absolute magnitudes and stellar 
,.1^ ■ masses. We consider both threshold and binned galaxy samples, with median B-band absolute magnitudes -21.6 < Mb - 5 log(/i) < -19.5 and 
^ I median stellar masses 9.8 < log(My,[/i"- Mg]) < 10.7. We assess the real-space clustering in the data from the projected correlation function, that 
we model as a power-law in the range 0.2 < r,,[/7"' Mpc ] < 20. Finally, we estimate the galaxy bias as a function of luminosity, stellar mass and 
redshift, assuming a flat A cold dark matter model to derive the dark matter 2PCF. 

Results. We provide the best-fit parameters of the power-law model assumed for the real-space 2PCF - the correlation length, ro, and the slope, 
^ [ 7 - as well as the linear bias parameter, as a function of the B-band absolute magnitude, stellar mass and redshift. We find a clear luminosity and 
^ stellar mass dependence of the galaxy clustering at all redshifts, with more luminous and massive galaxies showing a higher clustering amplitude, 
in agreement with previous studies. Our analysis provides new constraints on galaxy formation models that now have to match the observed 
redshift evolution of the clustering of VIPERS galaxies in a broad range of luminosity and stellar masses, between z = 0.5 and z= 1.1. 
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P\| , 1. Introduction through the mean occupation of galaxies in DM haloes, under 

' the assumption that all DM mass in the Universe is in the form 

' The Large Scale Structure (LSS) of the Universe contains fun- of haloes. 

^ : damental information to investigate the nature of the dark mat- i^^^^^ it has been observationally established in the 

ter (DM) and the origin of the accelerated expansion of the jo^^j Universe that the spatial distribution of galaxies 

Universe. However, since the LSS is traced by galaxies, the significandy depends on their luminosity, stellar mass, 

relevant information can only be extracted after specifying the ^^^^^^^ morphological and spectral types, and on the en- 

mapping between galaxies and DM. The way galaxies trace the yironment (see e.g. Davis & Geller 1976; Giovanelli et al. 

underlying density field. I.e. their biasing, can be understood ^ggg. Davis et al. 1988; Hamilton 1988; Einasto 1991; 

Maurogordato & Lachieze-Rey 1991; Loveday et al. 1995; 

Send ofprint requests to: F.Mamlli Benoist et al. 1996; Guzzo et al. 1997, 2000; Willmeretal. 

e-mail: federico.marulli3@unibo.it 1998; Brown et al. 2000, 2003; Norberg et al. 2001, 2002; 

* Based on observations collected at the European Southern Zehavi et al. 2002, 2005; Madgwick et al. 2003; Abbas & Sheth 

Observatory, Paranal, Chile under programs 182. A-0886 (LP) at the 2006; Zehavi et al. 2011; Budavari et al. 2003; Li et al. 

Very Large Telescope, and also based on observations obtained with 2006; Swanson et al. 2008; Loh et al. 2010; Ross et al. 201 1; 

MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, ^1. 2012; Christodoulou et al. 2012). The general trend 

at the Canada-France-Hawaii Telescope (CFHT), which is operated by ^^^^ luminous and massive galaxies with bulge-dominated 

the National Research Council (NRC) of Canada, the Institut National ui- jji . i.uf-. 

, „ . , . r 1 ^ XT • 1 J 1 Ti 1 1 morphologies and red colours cluster more strongly than faint, 

des Science de 1 Umvers of the Centre National de la Recherche . ^ • , , • , , • rr- i , 

Scientifique (CNRS) of France, and the University of Hawaii. This l^^^ massive, blue, spiral galaxies, an effect that becomes 

work is based in part on data products produced at TERAPIX and the particularly evident above the characteristic luminosity, L , and 

Canadian Astronomy Data Centre as part of the Canada-France-Hawaii stellar mass, M*, of the Schechter function. This is expected 

Telescope Legacy Survey, a collaborative project of NRC and CNRS. in the standard cosmological scenario of hierarchical growth of 

The VIPERS web site is http://vipers.inaf it/. cosmic Structures, where bright galaxies are hosted by massive, 
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rare and highly clustered DM haloes (White et al. 1987; Kaiser 
1987). 

Observations up to intermediate redshifts, z ~ I, i.e. half- 
way through the cosmic time from the microwave background 
to the present, show clustering trends similar to those ob- 
served in the local Universe, with the more luminous and mas- 
sive galaxies being more clustered than the faint, less massive 
ones (see e.g. Coil et al. 2006, 2008; McCracken et al. 2008; 
Meneux et al. 2006, 2008, 2009; Polio et al. 2006; Abbas et al. 
2010; Mostek et al. 2012). At higher redshifts, z > 1, observa- 
tional constraints are generally weaker, suffering from a com- 
bination of many types of selection biases, but the measure- 
ments still show a luminosity and stellar mass segregation in 
galaxy clustering (see e.g. Daddietal. 2003; Adelbergeret al. 
2005; Ouchi et al. 2005; Lee et al. 2006, 2009; Hildebrandt et al. 
2009; Wake et al. 2011; Lin et al. 2012), with strong evidence 
for a steepening of the correlation function at r < Ih^^ Mpc with 
increasing luminosity. 

It is possible to derive cosmological constraints from the 
spatial properties of galaxies using methods that do not explic- 
itly depend on galaxy bias, e.g. by exploiting baryonic acous- 
tic oscillations (see e.g. Eisenslein el al. 2005; Cole et al. 2005; 
Percival et al. 2007, 2010) or geometric clustering distortions 
(see e.g. Alcock & Paczynski 1979; Simpson & Peacock 2010; 
Marulli et al. 2012), or by combining the one- and two-point 
moments of the smoothed galaxy density distribution (Bel et al. 
2013). However, to extract the full information encoded in the 
shape of the galaxy power spectrum and in its redshift-space dis- 
tortions, the galaxy bias must be accounted for to some extent. 
Since the connection between baryons and DM can be signifi- 
cantly affected by non-linear phenomena, such as mergers, dy- 
namical friction, cooling, feedback, etc., semi-analytic models 
or full hydrodynamic simulations provide useful tools to prop- 
erly investigate how galaxies trace the DM distribution. Once a 
mathematical description of the galaxy bias has been obtained, 
cosmological constraints can be derived by marginalising over 
the uncertain parameters describing the bias. Alternatively, ob- 
servational measurements of the galaxy bias, derived in a given 
cosmological framework, can be used to constrain the parame- 
ters of galaxy formation models. 

In this paper, we focus on the dependence of galaxy cluster- 
ing on luminosity and stellar mass, providing new observational 
constraints at redshift 0.5 < z < 1.1, from the first data release 
of the ongoing VIMOS Public Extragalactic Redshift Survey 
(VIPERS, Guzzo et al. 2013). The comparison between VIPERS 
clustering measurements with similar ones at lower and higher 
redshifts can constrain the cosmic evolution of the relationship 
between DM and galaxy properties, and therefore between grav- 
ity and cosmology on one side and processes associated with 
baryonic physics on the other side. The large number of objects 
with accurate spectroscopic redshifts and the large volume cov- 
ered by the VIPERS survey give us the opportunity to measure 
the clustering of galaxies over a wide range of luminosities and 
stellar masses, up to the cosmic epoch when the Universe was 
still decelerating. In particular, we provide constraints on the 
projected correlation function, that we model as a single power- 
law. This represents a reasonable, though not optimal, approx- 
imation at rp < 20/1"' Mpc, where is the direction perpen- 
dicular to the line-of-sight, and facilitates the comparison with 
previous studies. A more detailed interpretation of these data us- 
ing Halo Occupation Distribution (HOD) modelling is deferred 
to a future work (see also de la Torre et al. 2013). 

The plan of this paper is as follows. In Section 2 we give a 
general overview of the VIPERS data, describing in detail the 



luminosity and stellar mass selections adopted in this work. We 
investigate the spatial distribution of VIPERS galaxies measur- 
ing the two-point correlation function (2PCF) in approximately 
volume-limited samples of different B-band absolute magnitude 
and stellar mass. The methodology of this analysis is described 
in Section 3, where we also explain how we account for the in- 
completeness affecting our VIPERS sample. More technical de- 
tails about the method used to correct for the proximity effect 
and stellar mass incompleteness are provided in Appendix A. 
The results of our analysis are presented in Section 4. Finally, in 
Section 5 we summarize the main results of this analysis. 

Throughout this paper, we assume a flat A cold dark mat- 
ter (CDM) background cosmology, with mass density parameter 
Q.M = 0.25, both to convert the measured redshifts into comov- 
ing distances and to estimate the DM correlation used to assess 
the galaxy bias. Changing the value of Qm in a reasonable range, 
given the current cosmological constraints, would introduce only 
small effects on our clustering measurements within their ob- 
servational uncertainties (see e.g. Marulli et al. 2012). Galaxy 
magnitudes are quoted in the AB system. The dependence of 
observed quantities on the Hubble parameter is indicated as a 
function of /? s //q/ 1 00 km s" ' Mpc" ' . 

2. The data 

2.1. The VIPERS survey 

We investigate the luminosity and stellar mass dependence of 
galaxy clustering using the VIPERS Public Data Release 1 
(PDR-1), that will be made publicly available in the fall of 2013. 
VIPERS is an ongoing Large Programme aimed at measuring 
redshifts for ~ 10^ galaxies at 0.5 < z < 1.2. One of the main 
goals of this survey is to characterize the galaxy spatial prop- 
erties and growth rate of the LSS, measuring galaxy clustering 
and redshift-space distortions at an epoch when the Universe 
was about half its current age. Indeed, according to the standard 
ACDM scenario, the redshift range that VIPERS probes is previ- 
ous to the dark energy (DE) dominated era, crossing the cosmic 
epoch when the Universe started accelerating (see the bottom 
labels in Fig. 1). 

The galaxy target sample is selected from the Canada- 
France-Hawaii Telescope Legacy Survey Wide (CFHTLS-Wide) 
optical photometric catalogues (Mellier et al. 2008). VIPERS 
covers ~ 24 deg- on the sky, divided over two areas within the 
Wl and W4 CFHTLS fields. Galaxies are selected down to a 
flux limit of < 22.5, and a robust gri colour pre-selection 
is applied to effectively remove objects at z < 0.5. Coupled 
with an aggressive observing strategy (Scodeggio et al. 2009), 
this allows us to double the galaxy sampling rate in the redshift 
range of interest (~ 40%), with respect to a pure magnitude- 
limited sample. At the same time, the area and depth of the sur- 
vey result in a fairly large volume, 5x10^ h^^ Mpc^ , analo- 
gous to that of the 2dFGRS at z ~ 0.1. Such a combination of 
sampling and depth is unique among current redshift surveys at 
z > 0.5. VIPERS is performed with the VIMOS multi-object 
spectrograph at the European Southern Observatory Very Large 
Telescope (Le Fevre et al. 2002, 2003), at moderate resolution 
{R = 210), using the LR Red grism, providing a wavelength 
coverage of 5500-9500A and a typical radial velocity error of 
140 km sec ' . The full VIPERS area is covered through a mo- 
saic of 288 VIMOS pointings (192 in the Wl area and 96 in the 
W4 area). A complete description of the VIPERS survey, from 
the definition of the target sample to the actual spectra and red- 
shift measurements, is presented in the parallel survey descrip- 



2 



F. Marulli et al.: Clustering as a function of luminosity and stellar mass 



comoving distance [h ' Mpc] comoving distance [h ' Mpc] 

1500 2000 2500 1500 2000 2500 




0.4 0.6 0.8 1 1.2 0.4 0.6 0.8 1 1.2 

redshift redshift 



Fig. 1. Luminosity and stellar mass selection as a function of redshift. The black dots show the Wl and W4 VIPERS galaxies. Yellow 
lines surround the selected sub-samples. Solid lines show the 90% completeness levels for the whole galaxy population (cyan lines) 
and for galaxies of four different photometric types: red E/Sa (red lines), early spiral (magenta lines), late spiral (green lines) and 
irregular or starburst (blue lines). The B-band rest-frame magnitude limits in the left panel account for the average luminosity 
evolution of galaxies, as explained in the text. In the bottom axes, Zeq,A~ 0.44 is the redshift of DM-DE equality and Zacc~ 0.82 
is the redshift when the Universe started accelerating, according to ACDM predictions with the set of cosmological parameters 
adopted in this work. 



tion paper (Guzzo et al. 2013), while a discussion of the survey 
data reduction and management infrastructure can be found in 
Garilli et al. (2012). 

The PDR-1 galaxy catalogue includes 53,608 redshifts and 
corresponds to the reduced data frozen in the VIPERS database 
at the end of the 2011/2012 observing campaign. The 2PCF 
of PDR-1 galaxies and the first constraints on the growth rate 
from redshift-space distortions are presented in a parallel paper 
(de la Torre et al. 2013). A measurement of Qm using a novel 
clustering statistic is given in Bel et al. (2013). The galaxy stellar 
mass function is presented in Davidzon et al. (2013), while the 
luminosity function and the evolution of the colour bimodality 
is provided in Fritz et al. (2013). The PDR-1 galaxies are classi- 
fied with a Support Vectors Machine algorithm in Malek el al. 
(2013), while an early subset of the spectra has been ana- 
lyzed and classified through a Principal Component Analysis in 
Marchetti et al. (2013). Finally, the VIPERS redshift distribution 
has been used to measure the real-space galaxy power spectrum 
with colour-selected samples from the CFHTLS (Granett et al. 
2012). 

2.2. Luminosity and mass selection 

We consider galaxy samples selected in three redshift ranges: 
z 6 [0.5,0.7], z e [0.7,0.9] and z e [0.9, 1.1]. To investigate the 
luminosity and mass dependence of galaxy clustering, we use 
both threshold and binned sub-samples with different selections 
in B-band absolute magnitudes, Mb, and stellar masses, M*. The 
latter quantities have been estimated with a spectral energy dis- 
tribution (SED) fitting technique, using the HYPERZMASS code 
(Bolzonella et al. 2000, 2010). Absolute magnitudes are mea- 



sured using the B Buser filter in the AB system. More details 
about the SED fitting procedure are given in Davidzon et al. 
(2013). 

Fig. 1 shows the sub-samples selected in this work. The 
adopted limits for the different sub-samples are highlighted by 
the grey area and the yellow solid lines. The black dots show 
the values of Mb {left panel) and M* (right panel) for the Wl 
and W4 VIPERS galaxies, as a function of redshift. The B-band 
magnitude limits account for the average luminosity evolution of 
galaxies. Following Meneux et al. (2009), we model the average 
redshift evolution of galaxy magnitudes as Mb(z) - Mb(0) + z 
(llbert el al. 2005; Zucca el al. 2009). Instead, for the stellar 
mass-selected sub-samples, we adopt the more conservative and 
widely used flat thresholds, as the redshift evolution of M* is 
negligible up toz ~ 1 (Pozzelli el al. 2007, 2010; Davidzon el al. 
2013). The solid coloured lines represent the 90% complete- 
ness levels for the whole galaxy population (cyan lines) and for 
galaxies of four different photometric types: red E/Sa (red lines), 
early spiral (magenta lines), late spiral (green lines) and irregu- 
lar or starburst (blue lines) (see Davidzon el al. 2013, for more 
details). We do not use these completeness limits to investigate 
the dependence of galaxy clustering on spectral type, that we de- 
fer to a future work. These lines are displayed here to show that 
our magnitude-selected samples are all volume-limited irrespec- 
tive of their colours, with the exception of the last redshift bin, 
where, due to our / band selection, we start losing the reddest 
galaxy population. However, the mass-selected samples are af- 
fected by the incompleteness of red galaxies. We will come back 
to this point in Section 3.4 and Appendix A. 

Finally, in the bottom axes of Fig. 1, we mark Zeq.A. the 
redshift of DM-DE equality, and Zacc^ the redshift when the 
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Universe started accelerating, according to ACDM predictions 
with the set of cosmological parameters adopted in this work. 
Clearly, the three redshift ranges considered for the present anal- 
ysis probe different cosmic epochs, crossing the transition from 
deceleration to acceleration. 

3. The two-point correlation function 

In this section, we describe the methodology used to measure 
the galaxy 2PCF. In particular, we discuss our error estimates 
and sampling rate correction. Finally, we describe the method 
used to correct for mass incompleteness. 

3.1. Clustering measurements 

We quantify the spatial properties of the VIPERS galaxies us- 
ing the 2PCF, ^(r). This function is implicitly defined as dP - 
n^[l + ^(r)]dVidV2, where dP is the probability of finding a pair 
of galaxies in the comoving volumes dVi and dV2, separated by 
a comoving distance r, and n is the galaxy average number den- 
sity. In the following, we will refer to real-space and redshift- 
space three-dimensional separations using the vectors r and s, 
respectively. 

To measure the function ^(i), we make use of the standard 
Landy & Szalay (1993) estimator: 

GG(s)-2GR(s) + RR(s) 

^^'^^ m^) ' 

where GG(s), GR{s) and RR{s) are the normalized galaxy- 
galaxy, galaxy-random and random-random pairs, with spatial 
separation in the range [s - As/2, s + Ai/2], respectively. 

The random samples are generated in the Right Ascension- 
Declination plane, adopting the same geometric mask of the real 
data. We set Nr - 3QNgai, where Nr is the number of random 
objects and Ngai is the number of galaxies in each sample. The 
latter values are reported in the fifth columns of Tables (1- 
4). The redshifts of the random objects are extracted from the 
smoothed radial distribution of the W1H-W4 VIPERS galaxies, 
estimated in each luminosity and stellar mass sub-sample. Fig. 2 
shows the whole comoving radial distances distribution of the 
Wl (green histograms) and W4 (yellow histograms) VIPERS 
galaxies. To smooth the total Wl -I-W4 radial distribution, we use 
the non-parametric Vmax method (see de la Torre et al. 2013, for 
more details). The blue solid line in Fig. 2 shows the result of the 
smoothing. A similar result can be obtained by filtering the dis- 
tribution with a Gaussian kernel of cr = 150/?"' Mpc, as shown 
by the red line. The advantage of the Vmax method is that it does 
not require us to make assumptions about the kernel. We have 
verified that all the results presented in this paper do not vary 
significantly when changing the smoothing scheme adopted. 
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Fig. 2. Comoving radial distance distribution of the Wl and 
W4 VIPERS galaxy samples, shown with the green and yel- 
low histograms, respectively. The red and blue lines show 
the total W1-I-W4 smoothed distribution, obtained by filter- 
ing the observed distribution with a Gaussian kernel of cr = 
150/i ' Mpc and using the y^ax method, respectively. 

below the resolution limit (de la Torre & Peacock 2012) and 
thus to be able to simulate the faintest and less massive galax- 
ies observed in VIPERS. A detailed description of the method 
used to construct our mock galaxy catalogues can be found 
in de la Torre et al. (2013). The covariance error matrix is es- 
timated from the dispersion among the mock catalogues; 

(2) 

k=l 

where ^, is the mean value of the 2PCF, measured in each 
mock catalogue. The errors on the 2PCF are then obtained from 
the square root of the diagonal values of the covariance matrix, 
cr, = VC/;, that include both the Poissonian noise and the un- 
certainties due to the size of the volume explored, i.e. the sam- 
ple variance. The luminosity and stellar mass dependence of the 
2PCF errors is estimated measuring the covariance matrix in 
mock sub-samples with the same luminosity and stellar mass 
selection used for the real data. 



3.2. Error estimate 

The errors on the clustering measurements are estimated using 
a set of mock galaxy catalogues. For the luminosity-selected 
samples, we use data-calibrated HOD mocks, while for the 
stellar mass-selected samples we consider galaxy mocks built 
from the stellar-to-halo mass relation (SHMR) of Moster et al. 
(2013). We use 26 independent HOD and SHMR mock cata- 
logues, both constructed by assigning galaxies to DM haloes 
of the MultiDark simulation, a large N-body run in the ACDM 
framework (Prada el al. 2012). A novel technique has been used 
to repopulate the MultiDark simulation with haloes of mass 



3.3. Weigliting tlie galaxy pairs 

To properly account for the effects of samphng rate and colour 
selection, each galaxy has been weighted according to its sky 
position and redshift, using the following weight function; 

W(Q,Z) - WCSR ■ WtsR ■ WSSR 

= CSR-\z)TSR-\Q)SSR-\Q). (3) 

All of the weights in Eq. (3) are derived directly from the data. 
The Colour Sampling Rate, CSR, represents the completeness 
due to the colour selection. It mainly depends on redshift and 
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can be modelled as CSR{z) = 0.5[1 - erf(7.405 - 17.465 ■ z)l 
where erf is the error function eTi{x) - 2/-\/7r^ e\p(-t^)dt 
(Guzzo et al. 2013; Coupon et al. 2013). The Target Sampling 
Rate, TSR, is defined as the fraction of galaxies in the pho- 
tometric catalogue that have been spectroscopically targeted, 
TSR - NgaiiO < Zflag ^ 9.5)/Marget, whcrc Zflag Is the quahty 
flag of the spectroscopic redshift measurement (see Guzzo et al. 
2013, for more details on Zflag)- In addition, we have a 3.2% stel- 
lar contamination in the target sample, but the effect on the TSR 
is negligible. Finally, the Spectroscopic Success Rate, SSR, is 
defined as the fraction of the spectroscopically targeted galaxies 
for which a secure identification (i.e. high quality flag) has been 
obtained, SSR = Ng^i(2 < za^g < 9.5)/A^gai(0 < Zflag < 9.5). Both 
the TSR and SSR depend on the VIMOS quadrant, Q, in which 
the galaxy has been observed. All the other possible dependen- 
cies, e.g. on redshift and magnitude, are comparatively smaller 
and can be neglected, as discussed in detail in de la Torre et al. 
(2013). 



expected to introduce a scale-dependent clustering suppression, 
up to ~50% at scales < 1 /i"' Mpc , for the redshift-space 2PCF. 

This correction assumes that different types of mock galax- 
ies have clustering properties that match those of the real ob- 
jects. However, this is only partially true. It is known that the 
model considered here overpredicts the abundance of faint and 
red galaxies. In addition, while the clustering of blue galaxies 
is well reproduced, the mock red galaxies appear more clus- 
tered than in real data (de la Torre et al. 2011; Cucciati et al. 
^. , _). For the above reasons, such a model-dependent correc- 
tion should be applied with caution. In the following, we will 
show our clustering measurements without correcting for stellar 
mass incompleteness. Nevertheless, in Tables 3-4, we provide 
our results obtained both with and without this correction, to di- 
rectly show what is the expected effect according to the semi- 
analytic model considered. 

More details about these two corrections for the proximity 
effect and mass incompleteness are provided in Appendix A. 



3.4. Proximity effects and mass incompleteness 

The VIMOS multi-object spectrograph allows spectra for ~400 
galaxies to be obtained in a single pointing (Bottini et al. 2005); 
however, it is not possible to target galaxies that are close to 
each other along the dispersion direction. Thus, the highest den- 
sity regions are under-sampled, resulting in a systematic under- 
estimation of the correlation function, especially on small scales 
(Polio et al. 2005; de la Torre et al. 2013). To estimate the im- 
pact of the proximity effect, we measure the mean ratio between 
the redshift-space 2PCF in two sets of 26 mock catalogues with 
and without the slit mask target selection algorithm applied. A 
similar test has been applied by Coil et al. (2008). The mock cat- 
alogues used here are constructed with a semi-analytic model 
(De Lucia & Blaizot 2007) on top of the Millennium simulation 
(Springel et al. 2005). To obtain a smoothed version of this ratio, 
we model it with an error function (see Fig. A. 2 in Appendix A). 
Then, we multiply the measured redshift-space 2PCF by this 
smoothed ratio. As we find this correction to be almost indepen- 
dent on luminosity and stellar mass, we use the values obtained 
in our largest galaxy sample to correct all of our clustering mea- 
surements. At z < 0.9, the under-estimation of the 2PCF due 
to the proximity effect is ~ 10% at scales r > 1 /i"' Mpc, and 
increases slightly at smaller separations. 

As shown in the right panel of Fig. 1, our mass-selected 
galaxy samples can be partially incomplete due to the limiting 
Iab < 22.5 flux cut of VIPERS, missing galaxies with high mass- 
to-light ratio. To correct for this effect, we would have to know 
what are the clustering properties of the galaxies that are missing 
in each of our samples. This can only be inferred using galaxy 
formation models both complete in mass and realistically de- 
scribing the formation and evolution of faint galaxies. 

Following what has already been done in previous studies 
(see e.g. Meneux et al. 2008, 2009), we implement a correc- 
tion method based on a similar strategy as described above for 
the proximity effect. Specifically, the correction is derived by 
measuring the ratio between the mean redshift-space 2PCF of 
two sets of mock catalogues, one complete in mass (down to 
~ IO^Mq) and the other one magnitude-limited at Iab < 22.5. 
According to the semi-analytic model considered here, the high 
mass-to-light ratio galaxies that are missed in our samples are 
faint red objects, that are predicted to be more clustered than the 
average population. Hence, the stellar mass incompleteness is 



3.5. Projected and real-space correlations and galaxy bias 

To correct for redshift-space distortions due to galaxy pecu- 
liar velocities, and derive the 2PCF in real-space, ^(r), we first 
measure the clustering as a function of the two distances per- 
pendicular and parallel to the line-of-sight, ^(rp,7T), using the 
Landy & Szalay (1993) estimator Then, we integrate along the 
line of sight, tt, as follows; 



Wp(rp) = 21 dn' ^{Vp, n') . 
Jo 



(4) 



The real-space 2PCF can then be estimated by deprojecting the 
function Wp{rp) (Davis & Peebles 1983; Saunders et al. 1992), 
as follows: 



^(0 



, dwp{r'p)/drp 



4^ 



(5) 



To facilitate the comparison with previous studies, we do not use 
Eq. 5, but model the function ^(r) as a single power-law: 



(6) 



With the above assumption, Eq. (4) can be solved analytically: 

(7) 



Wpirp) = rp\ — \ 



where F is the Euler's Gamma function. Hence, a power-law fit 
to Wpivp) provides ro and y for the real-space correlation func- 
tion ^(r). Finally, the scale-dependent galaxy bias can be esti- 
mated as: 



l^m(r)' 



(8) 



where ^m(r) is the matter 2PCF, that we derive using the software 
CAMB (Lewis & Bridle 2002), which exploits the HALOFIT rou- 
tine (Smith et al. 2003) to reproduce the non-linear evolution of 
the matter correlation function ' . 



Alternatively, the bias function can be also estimated directly from 
the projected correlation function, as a function of the direction perpen- 
dicular to the line-of-sight, i.e. birp) = ^Wp(rp)/w^(rp). 
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Fig. 3. 2PCF of VIPERS galaxies as a function of B-band absolute magnitude (upper panels) and stellar mass (bottom panels). The 
error bars are the square root of the diagonal values of the co variance matrix given by Eq. 2. Absolute magnitudes are given in the 
AB system. 



4. Results 

In this section, we present our results on the luminosity and stel- 
lar mass dependence of galaxy clustering and galaxy bias, and 
compare our results with those of previous studies. 

4.1. Redshift-space and projected 2PCF 

Fig. 3 shows the redshift-space 2PCF of VIPERS galaxies, ^(s), 
as a function of B-band absolute magnitude (upper panels) 
and stellar mass (bottom panels), in the three redshift ranges 
z G [0.5,0.7], z G [0.7,0.9] and z g [0.9, I.I]. The function 
^(s) has been estimated using Eq. (1), while the error bars are 
the square root of the diagonal values of the covariance matrix, 
given by Eq. (2). 

To derive the real-space 2PCF and directly compare the 
VIPERS clustering with previous data, we estimate the pro- 
jected correlation function using Eq. 4. We integrate ^(rp,n:)up 
to TTmax = 30 /i ' Mpc . As we have verified, this choice repre- 



sents a convenient compromise between robustness and the need 
to exclude noisy bins at large separations. The projected 2PCF of 
VIPERS galaxies is shown in Fig. 4. The upper and lower panels 
show, respectively, the luminosity and stellar mass dependence 
of Wp(rp), as indicated by the labels. The solid lines show the 
best-fit power-law models, given by Eq. 7, that we obtain fitting 
the correlation functions in the interval 0.2 < rp[/! ' Mpc ] < 20. 
The results presented in this work are robust with respect to the 
scale range used for the fit, except for the most luminous and 
massive galaxy samples. We can see in Fig. 4 that the power- 
law model does not describe well the small-scale behaviour of 
Wp(rp) for the most massive galaxies. This steepening in the am- 
plitude of the correlation function on scales below l/i ' Mpc can 
be interpreted in the framework of the halo model as due to the 
larger fraction of satellite galaxies in the massive haloes, creat- 
ing a prominent one-halo term in the correlation function. 

We find a clear luminosity and stellar mass segregation in 
galaxy clustering at all redshifts considered. This result is in 
agreement with previous findings, as discussed in detail in the 
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Fig. 4. Projected correlation function of VIPERS galaxies as a function of B-band absolute magnitude (upper panels) and stellar 
mass (bottom panels). The solid lines show the best-fit power-law models, obtained by fitting the correlation functions in the interval 
0.2 < rp[h-^ Mpc] < 20. 



next section, and provides new strong constraints on galaxy for- 
mation models in the intermediate redshift Universe. The clus- 
tering measurements presented in Fig. 3 and 4 refer to threshold 
galaxy samples, as indicated by the labels. Our measurements on 
the galaxy clustering in both threshold and binned galaxy sam- 
ples are given in Tables 1-4. Specifically, we provide the best- 
fit parameters of the power-law model: the correlation length ro 
(sixth columns), and the slope y (seventh columns). The errors 
on the best-fit parameters are derived from the scatter among the 
ro and y values measured in the 26 independent mock VIPERS 
mocks. 

4.2. Comparison to previous studies 

In Fig. 5, we compare our VIPERS measurements with many 
previous data, both at low and at high redshift. The upper pan- 
els show the best-fit values of ro and y as a function of the B- 
band magnitude. Following Polio el al. (2006), we normalize the 
median magnitudes of each galaxy sample to the corresponding 



values of the characteristic magnitude of the Schechter luminos- 
ity function, M* , obtained by fitting the values of M* (z) given 
by llbert et al. (2005). With this method, we can properly com- 
pare samples at different redshifts, taking into account the mean 
brightening of galaxies due to their evolution. 

Blue squares, red diamonds and green circles refer to our 
VIPERS constraints estimated at z 6 [0.5,0.7], z G [0.7,0.9] 
and z e [0.9, 1.1], respectively. A clear luminosity evolution of 
the clustering length, ro, is evident at all redshifts. We can see 
that ro increases with time between z ~ I and z ~ 0.6. Moreover, 
there is some evidence that ro grows faster for faint galaxies with 
respect to the bright ones. 

The black and grey symbols of Fig. 5 show previous mea- 
surements of the same quantity at different redshifts. Solid sym- 
bols refer to data taken in redshift ranges similar to ours, while 
open symbols are for data at lower or higher redshifts. For th e lo- 
cal Universe, we consider the data from Zehavi et al. (201 1) and 
Norberg el al. (2002). Zehavi el al. (2011) measured the lumi- 
nosity dependence of galaxy clustering using the Sloan Digital 
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Fig. 5. Best-fit values of the coirelation length, ro, {left panels) and slope, y, {right panels), as a function of B-band absolute 
magnitude {upper panels) and stellar mass {bottom panels). VIPERS costraints, obtained in threshold samples at z € [0.5,0.7], 
z 6 [0.7, 0.9] and z e [0.9, 1.1] are shown, respectively, by blue squares, red diamonds and green circles, as indicated by the labels. 
The remaining grey symbols refer to measurements based on data from other surveys, reported here for comparison: upper panels 
- DSS-DR7 data from Zehavi et al. (201 1), 2dfGRS data from Norberg el al. (2002), VVDS-Deep data from Polio et al. (2006) and 
DEEP2 data from Coil et al. (2006); lower panels - SDSS-DR2 data'^from Li et al. (2006), VVDS-Deep data from Meneux et al. 
(2008), NMBS data from Wake et al. (2011), DEEP2 data from Mostek et al. (2012) and BzK data from Lin et al. (2012). 



Sky Survey Seventh Data Release (SDSS-DR7), with volume- 
limited galaxy samples up to z = 0.25. Their results are shown 
by black open diamonds. Norberg et al. (2002) used the 2dF 
Galaxy Redshift Survey (2dFGRS) to measure the galaxy clus- 
tering up to z - 0.15. Dark grey open dots are used in Fig. 5 
to display their measurements. Grey open rotated triangles and 
sohd squares are the VIMOS-VLT Deep Survey (VVDS) clus- 



tering measurements at 0.1 < z < 0.5 and 0.5 < z < 1.2, 
respectively (Polio et al. 2006). Finally, the grey solid triangles 
show the DEEP2 Galaxy Redshift Survey results at z ~ 1 from 
Coil et al. (2006). The latter two clustering measurements are at 
redshifts similar to VIPERS. 

The bottom panels of Fig. 5 show the stellar mass de- 
pendence of ro and y. We compare our results to data 
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from SDSS (Li et al. 2006), VVDS-Deep (Meneux et al. 2008), 
NMBS (Wakeetal. 2011), DEEP2 (Mostek et al. 2012) and 
BzK (Lin et al. 2012). The best-fit clustering lengths from 
SDSS-DR2 data in the local Universe (Li et al. 2006) have been 
fitted by Mostek et al. (2012) and shown here by black open di- 
amonds. The measurements at intermediate redshifts similar to 
VIPERS are shown by grey squares, that represent the VVDS- 
Deep results obtained by Meneux et al. (2008) at z ~ 0.85, by 
grey triangles, that show the DEEP2 measurements at z ~ 1 
from Mostek et al. (2012), and by grey soUd rotated triangles, 
that show the clustering measurements by Wake et al. (2011), 
using NMBS photometric redshifts in the range 0.9 < z < 1.3, 
obtained with a fixed slope y = 1.6. All the remaining grey sym- 
bols refer to higher redshift data. Open triangles and pentagons 
refer to the Wake et al. (201 1) data at the two redshift ranges 
1.2 < z < 1.75 and 1.6 < z < 2.2, respectively. Finally, the dark 
grey circles show the ro values estimated from the angular clus- 
tering analysis of BzK-selected star-forming galaxies at z ~ 2, 
for a fixed slope of y = 1.8. 

The general picture that emerges from these previous stud- 
ies is that the clustering correlation length depends significantly 
on both luminosity and stellar mass, increasing monotonically 
from faint and low massive objects to luminous and massive 
ones, while the clustering slope remains fairly constant. It is 
interesting to note that for the intrinsically brightest samples, 
Mb > M*g, the clustering amplitude ro does not evolve signif- 
icantly from z ~ 1 measurements of VVDS and DEEP2 to the 
local references of 2dFGRS and SDSS. The same behaviour is 
seen in VIPERS: we find little evolution in the most luminous 
bin. The VIPERS data allow us to track ro with precision over 
a range of luminosities from z = 0.5 - 1.1, and we find that 
the evolution becomes more significant at fainter luminosities. 
A possible interpretation of this result is that the most luminous 



galaxies at z ~ 1 were already almost as strongly clustered as 
they are now. With the three redshift slices of VIPERS we can 
observe the slow accumulation of these structures from z ~ 1 . 1 
toz~0.5. 

On the other hand, when it comes to dependence of cluster- 
ing on stellar mass, not all previous measurements are consistent. 
In particular, as it can be noted in the bottom left panel of Fig. 5, 
there is a big scatter between previous measurements, even when 
they probe similar redshift ranges. Clustering lengths at z ~ 1 
found in the VVDS-Deep and DEEP2 surveys are significantly 
lower not only than the local measurements from the SDSS- 
DR7, but also than the values found by NMBS at 0.9 < z < 2.2 
and by the BzK at z ~ 2. Our estimates of ro are higher than the 
VVDS-Deep and DEEP2 ones, in the stellar mass range consid- 
ered, while they fairly agree with all the other data. 

The differences are likely to arise from a combination of 
systematic effects. The surveys considered in Fig. 5 are char- 
acterised by different properties, e.g. the methods for mass esti- 
mation, the cosmological parameters adopted to assess the dis- 
tances and the selection criteria. The latter certainly introduce 
the most relevant effects, as different surveys can be affected 
by stellar mass incompleteness to varying degrees. In particu- 
lar, the VIPERS mass-selected samples contain on average much 
brighter galaxies than VVDS and DEEP2, due to VIPERS higher 
flux selection cut. 

Moreover, as the power-law model does not describe per- 
fectly the data, measurements obtained fitting over different 
ranges could produce different values of the fit parameters. This 
effect is expecially significant for the most massive bins, where 
the shape of Wp(rp) is the furthest from a power-law and, thus, 
most sensitive to the range, binning and choice of points used 
for fitting. Finally, the impact of sample variance has been esti- 
mated with different techniques, and it might be underestimated 
in some cases. 

The VIPERS clustering slope is almost constant, considering 
the estimated errors, around y ~ 1.7-1.8, with only tentative 
evidences of a slight luminosity and redshift evolution. This is in 
overall agreement with all previous data, except for the relation 
measured in VVDS-Deep. Using the data from the VVDS-Deep 
Survey, Polio et al. (2006) found a significant steepening of the 
correlation slope from y ~ 1.6 at < Mb >= 19.6, to y ~ 2.4 at 
< Mb >= -21.3, fitting Wp{rp) for separation < 10 /j"' Mpc in 
their high-redshift range 0.5 < z < 1.2 (see the grey squares in 
the upper right panel of Fig. 5). 

To directly compare with VVDS data, we measure the 
VIPERS 2PCF in the same redshift range used by Polio et al. 
(2005), 0.5 < z < 1.2. In this way, any possible effect due 
to the intrinsic evolution in the VVDS broad redshift slice is 
present also in our measurements. We consider four volume- 
limited samples, with the following B-band magnitude thresh- 
olds Mb(z = 1.2)-51og(/j) < -21,-21.25,-21.5,-21.75. The 
best-fit values of y are shown in Fig. 6, compared to the VVDS- 
Deep measurements in volume-limited samples (grey squares). 
To investigate the scale dependence of the clustering slope, we 
compare the y values obtained fitting Wp ('"/>) in the interval 
0.2 < rp[/i 'Mpc] < 20 (green points), and considering only 
small scales, rp < 1 /i ' Mpc (red points). The proximity effect 
caused by the one-pass strategy used in VIPERS is corrected 
with the method described in Section 3.3 and Appendix A. We 
are confident that the correction recovers the shape of the cor- 
relation function on small scales to within 5% (de la Torre et al. 
2013). 
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Fig. 7. Bias of VIPERS galaxies, averaged over the range [5-10] /i ' Mpc , as a function of B-band absolute magnitude (left panel) 
and stellar mass (right panel). The coloured symbols used for VIPERS measurements are the same as in Fig. 5. The black solid line 
in the left panel shows the fitted bias function by Zehavi et al. (201 1). 



When the 2PCF is fit on scales up to 20/1 ' Mpc , y appears 
almost constant at the value ~ 1.7, as in Fig. 6, while if we 
fit only the small scales, y results to be larger and slightly in- 
creasing, from y ~ 1.8 at (Mb - M*) ~ -0.3, to y ~ 2 at 
<Mb - M;) ~ -0.9. The VIPERS clustering slope is still lower 
than the VVDS one at high luminosities, though the mismatch is 
not large (below 2cr significance). 

In conclusion, we find that the steeping of the slope in VVDS 
can be partly attributed to the range of scales used in the fit. 
However, the root cause of the discrepancy with VIPERS is 
probably due to sample variance in VVDS. 

4.3. The galaxy bias 

In this subsection, we present our estimates of the large-scale 
bias of VIPERS galaxies. We infer the scale-dependent bias, 
b{r), using Eq. 8. Then, we average b(r) in the range 5 < 
r[h^^ Mpc ] < 20, where the galaxy bias is fairly constant. We 
assume a flat ACDM model with the matter density parameter 
Qm = 0.25, the baryon density parameter Qb = 0.045, the scalar 
spectral index «, - 1 and the power spectrum normalization 
erg - 0.8. A full HOD formalism would be necessary to properly 
describe the bias relation between galaxies and DM at different 
scales. However, at large scales, the bias obtained from the ratio 
of the galaxy correlation function to that of matter is consistent 
with that obtained using HOD modelling (see e.g. Zehavi et al. 
2011). 

In Fig. 7, we present the galaxy bias estimated in each thresh- 
old VIPERS sample, as a function of B-band absolute magni- 
tude (left panel) and stellar mass (right panel). The error bars 
show the scatter among the bias values measured in the mock 
catalogues. The values of the bias and its error in each VIPERS 
sample are reported in Tables 1-4. The galaxy biasing estimated 
here depends on the assumed cosmological model through both 



cTg and Om- The amplitude of galaxy bias scales linearly with 
CTg'. The latter is known with an error of ~2%, only just neg- 
ligible when compared to uncertainties in our bias parameter 
(~4%). Uncertainties on Q.m, which is constrained in the inter- 
val ~ 0.22 - 0.28, also have an impact. Varying its input value in 
that range induces a ~2% shift in the estimated bias. On the other 
hand, relative bias between different galaxy types (e.g. the ratio 
between any two bias values reported in the tables), is almost 
model-independent. 

Fig. 7 provides a compact view of the dependence of 
VIPERS clustering on luminosity and stellar mass. More lu- 
minous and massive galaxies appear more strongly biased than 
fainter and less massive ones, at all redshifts considered. At fixed 
luminosity and stellar mass, the bias function increases with red- 
shift, while the amplitude of the 2PCF decreases (see Fig. 5). 
In the left panel of Fig. 7, the black solid line shows the lumi- 
nosity dependence of the galaxy bias in the local Universe, es- 
timated by Zehavi et al. (201 1) using the SDSS-DR7 data (see 
also Norberg et al. 2001; Tegmarketal. 2004, for similar re- 
sults). We note that the shape of the local bias results are broadly 
consistent with our findings at intermediate redshifts, while the 
normalisation is lower as predicted by evolution models (see e.g. 
Tegmark & Peebles 1998; Blanton et al. 2000; Pillepich et al. 
2010). 

A more detailed investigation on the scale dependence of the 
VIPERS galaxy bias derived with the deprojection method (see 
e.g. Marulli et al. 2012) and with a full HOD approach is de- 
ferred to a future work, while the non-linearity of the biasing 
function is presented in Di Porto et al. (2013). 



4.4. Comparison to tlieoretical models 

To provide a first theoretical interpretation of our measurements 
in a cosmological context, we compare the VIPERS cluster- 
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Fig. 8. Redshift-space 2PCF of VIPERS galaxies compared to theoretical expectations, as a function of B-band magnitude (upper 
panels) and stellar mass {bottom panels). The coloured symbols used for VIPERS measurements are the same as in Fig. 3. The grey 
shaded areas show the 2PCF measured in two large galaxy mock catalogues, constructed with the Munich semi-analytic model. The 
2PCFs in different sub-samples are offset by 0.5 dex, starting from the lower luminosity and stellar mass samples, for visual clarity. 



ing presented in previous sections with predictions of a stan- 
dard galaxy formation model. The result is shown in Fig. 8, 
that compares the redshift-space 2PCF of VIPERS galaxies, 
measured in threshold samples of different luminosities (upper 
panels) and stellar masses (lower panels), with the theoretical 
predictions of the Munich semi-analytic model (Blaizot et al. 
2005; De Lucia & Blaizot 2007). The coloured symbols used for 
VIPERS measurements are the same as in Fig. 3, while the grey 
shaded areas show the 2PCF measured in two large galaxy mock 
catalogues. The correlations in different sub-samples are offset 
by 0.5 dex, starting from the lower luminosity and stellar mass 
samples, for visual clarity. The Munich semi-analytic model, al- 
ready introduced in Section 3.4 to calibrate the corrections for 
the proximity effect and stellar mass incompleteness, has been 
developed on the outputs of the Millennium Run, a large N- 
body DM simulation of the ACDM cosmological framework 
(Springel et al. 2005). 



As it can be seen, our measurements can be understood 
within the standard ACDM-Hsemi-analytic framework. Both the 
luminiosity and stellar mass dependence of the VIPERS clus- 
tering appear in overall good agreement with theoretical ex- 
pectations, though we note that the more luminous and mas- 
sive VIPERS galaxies are slightly more clustered than the semi- 
analytic ones, expecially at high redshift. However, we also re- 
mark that the value of erg - 0.9, assumed in the Millennium sim- 
ulation, is quite higher with respect to the most recent cosmolog- 
ical constraints (see e.g. de la Torre et al. 201 1, for more discus- 
sions). If the Millennium predictions were rescaled accordingly, 
the agreement would slightly worsen. A detailed comparison be- 
tween our measurements and theoretical predictions from differ- 
ent galaxy formation models is beyond the scope of the present 
paper, and it will be fully investigated in a future work. 
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5. Summary and conclusions 

In this work, we studied how galaxy clustering depends on lu- 
minosity and stellar mass, in the redshift range 0.5 < z < 1.1. 
To quantify the effect, we measured the redshift-space 2PCF of 
VIPERS galaxies in the first data release, PDR-1, and derived 
the real-space clustering and bias assuming a power-law form 
for the projected correlation function. Errors were estimated us- 
ing a set of mock HOD and SHMR catalogues mimicking the 
characteristics and the observational selections of the VIPERS 
dataset. 

The main results of this work can be summarised as follows: 

- We confirm with unprecedented precision at z = 0.5 - 1.1 
a monotonic increase of the clustering length ro, both as a 
function of B-band magnitude and stellar mass, in all the 
three redshift ranges considered. On the contrary, the clus- 
tering slope, J, appears quite constant when derived in the 
range 0.2 < rp[h-^ Mpc ] < 20. 

- The dependence of the clustering length ro on luminosity is 
steepest at high redshift. In terms of the evolution of ro with 
redshift, we find that there is only a weak trend for the high- 
est luminosity galaxies, while significant evolution in ro is 
found in the lower luminosity samples. 

- VIPERS measurements are broadly consistent with previous 
studies, while establishing a bridge from the local references 
to the high-redshift Universe. However, we find a difference 
at the 2-cr level with previous VVDS-Deep and DEEP2 mea- 
surements at z ~ 1; we attribute this disagreement mainly 
to different stellar mass incompleteness levels in these cata- 
logues with respect to VIPERS. 

- We provide an estimate of galaxy bias, averaged over the 
range [5 - 20] h^^ Mpc , as a function of luminosity, stellar 
mass and redshift, assuming a standard flat ACDM frame- 
work. The results are consistent with the predictions of hi- 
erarchical galaxy formation, although we will devote future 
work to testing evolution models in detail. 

The measurements of the 2PCF presented here represent a 
considerable advancement over past surveys at z ~ 1 . The large 
volume of VIPERS ensures that sample variance is subdomi- 
nant. Furthermore, we have modelled the measurement covari- 
ance utilising mock surveys that account for both astrophysical 
and observational effects in detail. The results constitute new, 
stringent constraints on the galaxy 2PCF over a broad range of 
luminosities and stellar mass to be used as references against 
which to test models of cosmology and galaxy formation at red- 
shifts 0.5 - 1.1. 
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Table 1. Properties of the selected VIPERS sub-samples in threshold luminosity bins. The best-fit values of the clustering normal- 
ization, ro, and slope, y, have been obtained fitting Wp{rp) in 0.2 < rp[h^^ Mpc] < 20. The errors are from the scatter among the 
HOD VIPERS mocks. The bias values assume a flat ACDM model with Qm = 0.25 and erg = 0.8. 
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range 


redshift 


Mb(z =1.1) - 51og(/i) 


Mb(z = 1.1) - 51og(ft) 




[Ir' Mpc ] 








[0.5,0.7] 


0.62 


< -19.5 


-19.87 


17473 


4.45 ± 0.21 


1.67 


±0.05 


1.19 ±0.07 


[0.5,0.7] 


0.62 


< -20.0 


-20.15 


12432 


4.81 ±0.22 


1.69 


±0.05 


1.27 ±0.07 


[0.5,0.7] 


0.62 


< -20.5 


-20.49 


7472 


5.22 ± 0.30 


1.73 


±0.06 


1.33 ±0.09 


[0.5,0.7] 


0.62 


< -21.0 


-20.86 


3599 


5.58 ±0.39 


1.81 


± 0.08 


1.37 ±0.10 


[0.5,0.7] 


0.62 


< -21.5 


-21.28 


1236 


6.24 ± 0.56 


1.79 


±0.13 


1.53 ±0.15 


[0.7,0.9] 


0.79 


< -20.0 


-20.41 


14442 


4.56 ± 0.22 


1.65 


± 0.04 


1.32 ±0.07 


[0.7,0.9] 


0.80 


< -20.5 


-20.68 


9469 


4.95 ± 0.23 


1.67 


±0.04 


1.41 ±0.07 


[0.7,0.9] 


0.80 


< -21.0 


-21.05 
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5.32 ± 0.29 


1.70 


±0.05 


1.48 ±0.08 


[0.7,0.9] 


0.80 


< -21.5 


-21.45 
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5.95 ± 0.38 


1.72 


± 0.09 


1.62 ±0.11 


[0.9,1.1] 


0.97 


< -20.5 


-21.00 


5018 


4.29 ±0.19 


1.63 


± 0.04 


1.37 ±0.06 


[0.9,1.1] 


0.98 


< -21.0 


-21.25 


3035 


5.08 ± 0.26 


1.64 


±0.05 


1.57 ±0.08 


[0.9,1.1] 


0.99 


< -21.5 


-21.65 


1123 


5.87 ± 0.43 


1.68 


± 0.08 


1.76 ±0.11 



Table 2. Properties of the selected VIPERS sub-samples in binned luminosity bins. Parameters and errors are as in Table 1. 
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redshift 


Mb{z= l.l)-51og(/7) 


Mb{z= l.l)-51og(/7) 


Mpc ] 








[0.5,0.7] 


0.62 


[-20.5,-19.5] 


-19.51 


10001 


3.92 ±0.19 


1.65 


± 0.05 


1.09 ±0.06 


[0.5,0.7] 


0.62 


[-21.0,-20.0] 


-19.95 


8833 


4.58 ±0.21 


1.66 


±0.05 


1.23 ± 0.07 


[0.5,0.7] 


0.62 


[-21.5,-20.5] 


-20.40 


6236 


5.00 ±0.29 


1.72 


±0.06 


1.29 ±0.09 


[0.5,0.7] 


0.62 


[-22.0,-21.0] 


-20.83 


3334 


5.35 ±0.36 


1.78 
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6.29 ± 0.59 


1.80 


±0.13 


1.53 ±0.14 


[0.7,0.9] 


0.78 


[-21.0,-20.0] 


-20.19 
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4.08 ±0.21 


1.61 


±0.05 


1.23 ±0.07 


[0.7,0.9] 


0.80 


[-21.5,-20.5] 


-20.59 
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4.71 ±0.24 


1.67 


±0.04 


1.35 ±0.07 


[0.7,0.9] 


0.80 


[-22.0,-21.0] 


-21.01 
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±0.05 
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0.96 
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-20.86 
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3.63 ±0.15 


1.64 


± 0.04 


1.19 ±0.05 


[0.9,1.1] 


0.98 


[-22.0,-21.0] 


-21.21 


2111 


4.90 ± 0.29 


1.65 


±0.04 


1.52 ±0.08 


[0.9,1.1] 


0.99 


[-22.5,-21.5] 


-21.64 


1091 


5.98 ± 0.43 


1.70 


± 0.08 


1.77 ±0.11 



Table 3. Properties of the selected VIPERS sub-samples in threshold stellar mass bins. The quantities r^, j'^ and bias'' are the 
clustering normalization, slope and linear bias corrected for the stellar mass incompleteness, respectively. The other parameters and 
errors are the same as in Table 1 . 
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range 


redshift 


[h-^ Mo] 


[h-^ Mo] 


[r'Mpc] 












[0.5,0.7] 


0.61 


> 9.0 


9.82 


17102 


4.55 ± 0.22 


1.69 


±0.05 


1.20 ±0.07 


4.99 


1.70 


1.30 


[0.5,0.7] 


0.62 


> 9.5 


10.11 


11567 


5.14 ±0.25 


1.76 


±0.05 


1.30 ±0.08 


5.70 


1.77 


1.42 


[0.5,0.7] 


0.62 


> 10.0 


10.35 


6880 


5.84 ±0.28 


1.83 


±0.05 


1.42 ±0.08 


6.01 


1.83 


1.46 
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> 10.5 
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1.86 


±0.06 


1.65 ±0.10 


6.97 


1.86 


1.66 
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> 9.0 


9.93 
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4.53 ± 0.25 
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[0.7,0.9] 
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>9.5 
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4.92 ± 0.26 


1.67 


±0.04 


1.40 ±0.08 


5.41 


1.69 


1.51 


[0.7,0.9] 


0.79 


> 10.0 


10.39 
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5.37 ± 0.29 


1.72 


±0.04 


1.48 ± 0.09 


5.87 


1.73 


1.59 


[0.7,0.9] 


0.79 


> 10.5 


10.67 


2498 


7.08 ±0.38 


1.73 


±0.05 


1.87±0.11 


7.22 


1.73 


1.90 


[0.9, 1.1] 


0.97 


> 9.5 


10.19 
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4.48 ±0.21 


1.60 


±0.04 


1.44 ±0.07 


4.51 


1.62 


1.44 


[0.9,1.1] 


0.97 


> 10.0 


10.46 


2857 


5.57 ±0.27 


1.62 


±0.04 


1.71 ±0.08 


5.78 


1.65 


1.74 


[0.9,1.1] 


0.97 


> 10.5 


10.71 
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6.52 ± 0.35 


1.74 


±0.05 


1.87 ±0.10 


6.78 


1.76 


1.93 
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Table 4. Properties of the selected VIPERS sub-samples in binned stellar mass bins. Parameters and errors are as in Table 3. 
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redshift 
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median stellar mass 
[h-- Mo] 
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7 
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/' 


bias*- 


[0.5,0.7] 
[0.5,0.7] 
[0.5,0.7] 


0.62 
0.62 
0.62 


[9.5, 10.5] 
[10.0, 11.0] 
[10.5,11.5] 


9.88 
10.19 
10.66 


7487 
2952 
2150 


5.47 
6.54 
6.93 


±0.31 
±0.32 
±0.36 


1.82 
1.87 
1.87 


±0.06 
±0.05 
±0.06 


1.34 ±0.09 
1.55 ±0.09 
1.64 ±0.10 


6.36 
6.73 
6.95 


1.84 
1.87 
1.87 


1.53 
1.60 
1.65 


[0.7,0.9] 
[0.7, 0.9] 
[0.7, 0.9] 


0.79 
0.79 
0.79 


[9.5, 10.5] 
[10.0, 11.0] 
[10.5, 11.5] 


9.92 
10.19 
10.67 


7410 
3353 
2497 


5.11 
6.31 
7.04 


±0.32 
±0.33 
±0.37 


1.69 
1.72 
1.73 


±0.06 
±0.05 
±0.05 


1.43 ±0.10 
1.70 ±0.10 
1.86±0.11 


6.25 
6.95 
7.21 


1.69 
1.73 
1.73 


1.70 
1.84 
1.90 


[0.9, 1.1] 
[0.9, 1.1] 


0.97 
0.97 


[9.5, 10.5] 
[10.0, 11.0] 


9.98 
10.25 


2942 
1620 


5.21 
6.23 


±0.35 
±0.30 


1.60 
1.77 


±0.06 
±0.05 


1.63 ±0.09 
1.78 ±0.09 


5.63 
6.52 


1.60 
1.78 


1.74 
1.85 
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Appendix A: Correcting for stellar mass 
incompleteness and proximity effect 

A possible method to correct the 2PCF measurements for stel- 
lar mass incompleteness is described in Section 3.4. We pro- 
vide here more details on this procedure and quantify the im- 
pact on our meaurements. The method consists in measuring the 
2PCF in two sets of mock galaxy catalogues, one set complete 
in stellar masses, while the other one with the same VIPERS 
lAB < 22.5 flux cut. We use mock catalgoues constructed with 
the Munich semi-analytic galaxy formation model (Blaizot et al. 
2005; De Lucia & Blaizot 2007), on top of the DM halo trees 
from the Millennium simulation (Springel et al. 2005). These 
catalogues are complete in stellar masses down to ~ 1O**M0. 
The red dots in Fig. A. 1 show the mean ratio between the 2PCF 
in flux-limited and complete mocks, averaged over the 10 inde- 
pendent catalogues, while the error bars are the rms scatter. The 
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Fig. A.2. As Fig. A.l, except for the redshift-space 2PCF mea- 
sured in mock catalogues with and without the slit mask target 
selection algorithm applied. 



Fig. A.l. Ratio between the redshift-space 2PCF measured in 
flux-limited and complete semi-analytic mocks. The red dots 
show the clustering ratio averaged over 10 independent mock 
catalogues, while the error bars are the rms scatter The black 
lines show the best-fit error function used to smooth these ratios. 



results refer to different redshift and stellar mass limits, as ex- 
plicitly indicated by the labels. According to the Munich model, 
the mass incompleteness introduces a clear scale-dependent re- 
duction of clustering, mainly significant at scales < 1 /i"' Mpc . 
The black lines, that show the best-fit error functions used to 
smooth these ratios, represent the correction that we actually use 
to obtain the values of , j'- and bias*- reported in Tables (3-4). 
As already stressed in Section 3.4, since the mock red and faint 
galaxies are significantly more clustered than the observed ones, 
this model-dependent correction should be taken with caution. 

As discussed in Section 3.3, due to our multi-object slit 
masks, it is not possible to observe galaxies that are too close 
to each other This proximity effect introduces a further scale- 
dependent suppression in the 2PCF that we have to correct. We 
adopt the same strategy described above to correct for mass in- 
completeness, i.e. we measure the ratio between the redshift- 
space 2PCF in mock catalogues with and without the slit mask 
target selection algorithm applied (also see Coil et al. 2008). As 
before, the mock galaxy catalogues used here are constructed 
with the Munich semi-analytic galaxy formation model. Fig. A. 2 
shows the results of this method. The red dots are the mean val- 
ues of this ratio as a function of luminosity, stellar mass and 
redshift, as reported in the labels. The error bars show the scat- 
ter between the mock samples. The black lines, as in Fig. A.l, 
show the error functions used to model these ratios. This correc- 
tion results to be almost independent on luminosity and stellar 
mass, similar to the results of Coil et al. (2006). Thus, we de- 
cide to use the correction calibrated in our largest galaxy sample 
to correct all our clustering measurements. We notice that the 
systematic clustering suppression due to the proximity effect be- 
comes less severe going to high redshift. More details can be 
found in de la Torre et al. (2013). 



16 



F. Marulli et al.: Clustering as a function of luminosity and stellar mass 



Dipartimento di Fisica e Astronomia - Universita di Bologna, viale 

Berti Pichat 6/2, 1-40127 Bologna, Italy 
^ INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, I- 

40127 Bologna, Italy 
^ INFN - Sezione di Bologna, viale Berti Pichat 6/2, 1-40I27 Bologna, 

Italy 

Dipartimento di Matematica e Fisica, Universita degli Studi Roma 
Tre, via della Vasca Navale 84, 00146 Roma, Italy 
5 INFN - Sezione di Roma Tre, via della Vasca Navale 84, 1-00146 
Roma, Italy 

* INAF - Osservatorio Astronomico di Roma, via Frascati 33, 1-00040 
Monte Porzio Catone (RM), Italy 

SUPA - Institute for Astronomy, University of Edinburgh, Royal 
Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 
INAF - Osservatorio Astronomico di Brera, Via Brera 28, 20122 
Milano, via E. Bianchi 46, 23807 Merate, Italy 
' Dipartimento di Fisica, Universita di Milano-Bicocca, P.zza della 
Scienza 3, 1-20126 Milano, Italy 

Astronomical Observatory of the Jagiellonian University, Orla 171, 
30-001 Cracow, Poland 
" National Centre for Nuclear Research, ul. Hoza 69, 00-681 
Warszawa, Poland 

INAF - Osservatorio Astronomico di Torino, 10025 Pino Torinese, 
Italy 

Aix Marseille Universite, CNRS, LAM (Laboratoire 
d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, 
France 

Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Highway, 
Kamuela, HI 96743, USA 

Centre de Physique Theorique, UMR 6207 CNRS -Universite de 
Provence, Case 907, F- 13288 Marseille, France 
INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via 
Bassini 15, 20133 Milano, Italy 
" Laboratoire Lagrange, UMR7293, Universite de Nice Sophia- 
Antipolis, CNRS, Observatoire de la Cote d'Azur, 06300 Nice, 
France 

'** Institute of Astronomy and Astrophysics, Academia Sinica, P.O. 

Box 23-141, Taipei 10617, Taiwan 
" INAF - Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 

34143 Trieste, Italy 

Institute of Physics, Jan Kochanowski University, ul. Swietokrzyska 
15, 25-406 Kielce, Poland 

Department of Particle and Astrophysical Science, Nagoya 
University, Furo-cho, Chikusa-ku, 464-8602 Nagoya, Japan 
^- Institute d'Astrophysique de Paris, UMR7095 CNRS, Universite 
Pierre et Marie Curie, 98 bis Boulevard Arago, 75014 Paris, France 
Universitatsstemwarte Miinchen, Ludwig-Maximillians Universitat, 
Scheinerstr. 1, D-81679 Miinchen, Gemany 

Max-Planck-Institut fiir Extraterrestrische Physik, D-84571 
Garching b. Miinchen, Germany 

Institute of Cosmology and Gravitation, Dennis Sciama Building, 
University of Portsmouth, Bumaby Road, Portsmouth, POl 3FX, 
UK 

INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Bologna, 
via Gobetti 101, 1-40129 Bologna, Italy 

INAF - Istituto di Radioastronomia, via Gobetti 101, 1-40129 
Bologna, Italy 

Universita degli Studi di Milano, via G. Celoria 16, 20130 Milano, 
Italy 



17 



